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Abstract 

We systematically investigate the parameter space of neutrino and charged lepton mass 
matrices for textures motivated by an extended quark-lepton complementarity. As the 
basic hypothesis, we postulate that all mixing angles in JJg and U u be either maximal or 
described by powers of a single small quantity e ~ 6q. All mass hierarchies are described 
by this e as well. In this study, we do not assume specific forms for Ug and U u , such as large 
mixing coming from the neutrino sector only. We perform a systematic scan of the 262 144 
generated mixing matrices for being compatible with current experimental data, and find 
a sample of 2 468 possibilities. We then analyze and classify the effective charged lepton 
and neutrino mass textures, where we especially focus on a subset of models getting under 
pressure for small 6 ri . In addition, we predict the mixing angle distributions from our sample 
of all valid textures, and study the robustness of this prediction. We also demonstrate how 
our procedure can be extended to predictions of the Dirac and Majorana phases in C/pmns- 
For instance, we find that CP conservation in neutrino oscillations is preferred, and we can 
impose a lower bound on the mixing matrix element for Qv/3(3 decay. 
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1 Introduction 



The Standard Model (SM) of elementary particle physics with three light Majorana neutri- 
nos contains 28 free parameters. Most of them, in total 22, describe the masses and mixings 
of the fermions (the remaining six parameters are the three SM gauge couplings, the vac- 
uum expectation value and quartic coupling of the Higgs, and the QCD 8 parameter). This 
large number of parameters, especially in the fermion sector, is generally considered as an 
unsatisfactory feature of the SM and one therefore seeks for models in which the number 
of parameters can be minimized. One possibility to reduce the number of parameters is to 
embed the SM into a Grand Unified Theory (GUT). 

By putting quarks and leptons into GUT multiplets, the masses and mixing angles in the 
quark and lepton sectors become related. From this point of view, it is thus reasonable to 
describe the observed hierarchical pattern of the masses and mixing angles of quarks and 
charged leptons [1] in terms of powers of a single small expansion parameter e. The expansion 
parameter e might, for example, represent a low-energy remnant of a flavor symmetry that 
has been broken at some high scale. In fact, the CKM mixing matrix Vckm [2,3] exhibits 
quark mixing angles of the orders 

\V w \~e, \V cb \~e\ \V ub \~e 3 , (1) 

where the quantity e is of the order of the Cabibbo angle 8q — 0.2. Similarly, for the same 
value e ~ 6*c, the mass ratios of the up- type quarks, down-type quarks, and the charged 
leptons can be approximated, e.g., by0 

m u : m c : m t = e 6 : e 4 : 1, m^ : m s : m^, = e 4 : e 2 : 1, m e : m M : m T = e 4 : e 2 : 1, (2) 

where mb/m t ~ e 2 , m T jm b ~ 1, and m t ~ 175 GeV. While the CKM angles and charged 
fermion masses are thus strongly hierarchical, there are striking differences in the neutrino 
sector. In the past few years, solar [4,5], atmospheric [6], reactor [7,8], and accelerator [9] 
neutrino oscillation experiments have established with increasing precision that among the 
leptonic mixing angles only the reactor angle #13 is small whereas the solar angle 812 and the 
atmospheric angle 823 are both large. Moreover, neutrino oscillation data tells us that the 
neutrinos have only a mild hierarchy. To be specific, expressing the neutrino mass ratios as 
in Eq. (T2]) in terms of powers of e, the neutrino mass spectrum can, e.g., be written as 

mi : m 2 : m 3 = e 2 : e : 1, m x : m 2 : m 3 = 1 : 1 : e, m 1 : m 2 : m 3 = 1:1:1, (3) 

where mi, m 2 , and m.3 denote the 1st, 2nd, and 3rd neutrino mass eigenvalue. In Eq. ([3]), the 
first equation corresponds to a normal hierarchical, the second to an inverted hierarchical, 
and the third to a degenerate neutrino mass spectrum^] In addition, we know from cosmolog- 
ical observations that the absolute neutrino mass scale is of the order ~ 10~ 2 . . . 10 _1 eV [10]. 
An attractive origin of the smallness of neutrino masses is provided by the seesaw mech- 
anism [11, 12]. For a summary of current values and errors for the neutrino oscillation 
parameters from a global analysis, see Table [H 

1 We take here a fit compatible with an SU(5) GUT. 

2 More precisely, we have in Eq. ([3]) for the inverse hierarchical case mi > mi and {mi — m\)lmi ~ e 2 . 
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Parameter 


Best-fit ±1 a 


la acc. 


2a range 


3cr range 


Projection ~ 2016 


Am 2 21 [10" 5 eV 2 ] 
lAmlil [10" 3 eV 2 ] 


7.9 ±0.3 


4% 


7.3-8.5 


7.1-8.9 




o c+0.20 
z - d -0.25 


10% 


2.1-3.0 


1.9-3.2 




sin 2 9 12 


o.so;^ 


9% 


0.26-0.36 


0.24 - 0.40 


4.6% (la) 


sin 2 9 23 


0.501S5? 


16% 


0.38-0.64 


0.34 - 0.68 


10% (la) 


sin 2 0i3 






< 0.025 


< 0.041 


< 0.0076 (3a) 



Table 1: Current best-fit values with la errors, relative accuracies at la, as well as 2a and 3tr allowed 
ranges of three- flavor neutrino oscillation parameters (from a combined analysis of all available data [13]). In 
addition, we show as the last column the projected improvement of the mixing angle errors on a time scale of 
about ten years from now. For this projection, we assume that the current best-fit values remain unchanged, 
in particular, that sin 2 29\3 will not be discovered. We use the experimental bounds sin 2 2#i3 < 0.03 (3a, 
year 2016 from Ref. [14] based on NO^A), ~10% precision of sin 2 623 (la, from Ref. [15], based on beam 
experiments), 4.6% precision of sin 2 9±2 (la, from Ref. [16] based on SPMIN/SADO+others with an exposure 
of lOGWktyr, which seems to be reasonable on that time scale for a KamLAND-scale detector). These 
values should be interpreted with care, because they depend on the experimental strategy. 

In spite of the qualitative differences between the quark and the lepton sector, there have 
recently been proposed interesting "quark-lepton complementarity" (QLC) relations [17-19] 
(for an early approach see Ref. [20]) that might point to quark-lepton unification. The QLC 
relations express that the solar angle 9 12 and the atmospheric angle #23 seem to be connected 
to the quark mixing angles by 

9 12 + #c ~ tt/4, 9 23 + 9 cb « vr/4, (4) 

where 9 cb = arcsin V^. A simple interpretation of how the QLC relations in Eq. might 
arise is to assume that the small and the maximal angles in this relation, 9q, 9 c b, and 
7r/4, describe the mixing of the left-handed charged leptons and the neutrinos in flavor 
basis. Consequently, the observed large leptonic mixing angles 9 12 and 6*23 can only arise 
as a result of taking the product of the charged lepton mixing matrix Ue and the neutrino 
mixing matrix U v in the PMNS [21,22] mixing matrix £7pmns = U\U U . This point of view 
is supported by the fact that in explicit models, 9q and 9 c b are usually given from the 
outset and also maximal leptonic mixing has meanwhile been obtained in many models (for 
a review of recent developments see, e.g., Ref. [23]). Therefore, if the QLC relations really 
arise from taking the product of U e and U v in [/pmnsj it is interesting to study models with 
mass matrices where, e.g., the value of 9 12 can be understood in terms of 9c and #23 as a 
result of combining Ue and U v into t/pMNS- 

There has been a considerable amount of work on QLC. For example, deviations from 
bimaximal neutrino mixing [24] that are introduced by 9q have been studied in Refs. [25- 
27]. Sum rules in the context of QLC using some assumptions on the mixing angles were 
presented in Refs. [28,29]. Parameterizations of £7pmns i n terms of 9q as an expansion 
parameter were given in Refs. [30-34]. For renormalization group effects and QLC see 
Ref. [35] and references therein. Model building realizations of QLC were discussed in 
Refs. [36,37] and for the special case 9c — 6*13 in Ref. [38]. 

In this paper, we study systematically 262 144 pairs of charged lepton and neutrino mass 
matrices that lead to the observed leptonic masses and mixing angles. Motivated by QLC, 
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our only assumption on the mixing angles in Ui and U v is that they can take any of the 
values 7r/4, e, e 2 ,0, where e is of the order of the Cabibbo angle e ~ 6q and represents the 
only small expansion parameter in our approach. Different from earlier studies of QLC, 
which assume very specific forms of Ug and U u [25-29], we allow maximal (~ 7r/4) and small 
(~ e, e 2 , . . . ) mixing angles to be generated in the charged lepton as well as in the neutrino 
sector. This means, in particular, that we do not assume that Ug or U v are necessarily of a 
bimaximal mixing form but obtain this, instead, as a special case. We consider all possible 
combinations of the mixing angles 7r/4, e, e 2 , . . . in Ut and U u , and select those pairs Ui and 
U u that give PMNS mixing angles #12, #13, and #23 in agreement with observations. For the 
detailed description of our method, see Sec. [2j We then unambiguously reconstruct directly 
the charged lepton and the neutrino mass matrices from the known charged lepton and 
neutrino mass spectra in Eqs. (j2J) and (j3J). The matrix pairs that we obtain in this way 
are then presented in the flavor basis with order of magnitude entries that are expressed 
as powers of e ~ 9q, which takes the role of a small expansion parameter for the matrices. 
Such approximate forms or representations of mass matrices expanded in terms of e will be 
called in the following "mass matrix textures" or simply "textures" . They are presented and 
discussed in Sec. |3j We analyze the predictions for the mixing angles from the set of mixing 
matrices consistent with current observations in Sec. HI and we study the robustness of these 
"predictions". Finally, we show how our framework can be applied to phase predictions in 
Sec. [5j such as for CP violation in neutrino oscillations or the mixing matrix element in 
neutrinoless double beta decay. 



2 Method 



In this section, we will introduce our method for a systematic scan of the textures that 
exhibit QLC. First, we set up the notation for the mass and mixing matrices. Next, we 
motivate and present our hypothesis for the charged lepton and neutrino mixing angles. 
We show how the resulting PMNS matrices are compared with observation using a selector 
and discuss how to construct from the selected mixing matrices the textures for the charged 
leptons and neutrinos. 



2.1 Mixing Formalism and Notation 



In this section, we describe the formalism and notation for the leptonic mixing parameters. 
Here, we follow closely Ref. [27]. A general unitary 3x3 matrix U un i tar y can always be 
written as 

Unitary = diag (e^ 1 , e^ 2 , )-U- diag (e iQl , e iQ2 , l) , (5a) 
where the phases tpi, <f 2 , <f3, and a> 2 , take their values in the interval [0, 2%] and 

/ C12C13 S19C13 si^e -1 ^ \ 



U 



-S12C23 - Ci 2 S 2 3Sl3e 

Ad 



iS 



S12C13 
C12C23 - S12S23S13& 



id 



S23C13 



\ S12S23 - Ci2C23Si 3 e 10 -C12S23 - Si2C23Si 3 e 1 ' 5 C23C13 / 



(5b) 
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is a CKM-like matrix in the standard parameterization with Sij = sin 9ij, Qj = cos 9ij, where 
9ij G {On, #23} he all in the first quadrant, i.e. §ij G [0, f ] , and 5 G [0, 27r]. The matrix 
Z7 is thus described by 3 mixing angles 9ij and one phase 5, i.e., it has 4 parameters. The 
matrix t/ un itary has five additional phases and contains therefore in total 9 parameters. 

Depending on whether neutrinos are Majorana or Dirac particles, the low-energy effective 
Lagrangian for lepton masses takes one of the forms 

C M = -(M^e) - l(Mf*%^ + h.c, £ D = -(M t ) ij e i e D j - (M™™)^ + h.c, (6) 

where and Ui are the left-handed charged leptons and neutrinos that are part of the SU(2) 1 
lepton doublets £j = (vi, e;) T , while ef are the right-handed (SU(2)l singlet) charged leptons, 
and i = 1,2,3 is the generation index. We have also extended the SM by adding to each 
generation i one right-handed SM sing let neutrino uf. In Eq. (jgj), M e , M^ aj , and M° irac , 
denote the charged lepton (Me), the Majorana neutrino (M;7 aJ ), and the Dirac neutrino mass 
matrix (M° irac ), respectively. The Dirac matrices Me and M° irac are general complex 3x3 
matrices, whereas M^ 1 ^ is a complex symmetric 3x3 matrix. The charged lepton mass 
matrix is diagonalized by a biunitary transformation 

Me = UeMf ag U'e j , (7a) 

where Ue and U' e are unitary matrices acting on the left-handed (Ue) and right-handed (U' £ ) 
charged leptons and ef, which span the rows and columns of Me, respectively. In Eq. (I7a|) . 
the matrix M diag is on the diagonal form M diag = diag(m e , m M , m T ). Using the freedom of 
re-phasing the charged lepton fields, we can, in what follows, assume that Ui is on a CKM- 
like form that is parameterized as in Eq. (15bj) . We then define the 4 mixing parameters 
of the left-handed charged leptons, the three mixing angles 9[ 2 , 9{ 3 , 9 23 , and the phase 5 e , 
by identifying in Eq. f!5bj) 9^ — > 9\- and 5 —>■ 5 e . The Majorana and Dirac neutrino mass 
matrices are given by 

M, Maj = U u Mf^Ul, M° irac = U„ M, diag U'J, (7b) 

where the unitary neutrino mixing matrix U u acts on the left-handed neutrinos z/j, while 
U' v acts on the right-handed neutrinos vf. In Eq. (17bft . M diag is a diagonal matrix M diag = 
diag(mi, rn.2, ^3). Since Ue has already been brought to a CKM-like form, we have no longer 
the same freedom to remove phases in U u and, thus, we find that the PMNS matrix is in 
general written as 

t/pMNS = U\U V = U\DU V K, (8a) 

where U v is a CKM-like matrix that is on the form as in Eq. (15bj) while D = diag(l, e lipi , e 1{p2 ) 
and K = diag^ 1 ^ 1 , e 1 "^ 2 , 1) are diagonal matrices with phases in the range <pi,<f2,<l>i,<j>2 G 
[0, 27r]. Note that we have already removed in £/pmns an unphysical overall phase. The 
CKM-like matrix U u in Eq. flHaj) contains four neutrino mixing parameters, the three neu- 
trino mixing angles 9 v vl ,9 v Vi ,9Y l?) , and a neutrino phase 8 V , which we define in the standard 
parameterization by identifying in Eq. (j5bjl the neutrino mixing angles as 9^ — > 9^ and the 

neutrino phase as 5 — > 5 U . The matrix in Eq. f|5a|) is written in terms of the six phases 
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8 , 8 U , tpi, if 2, 0i, and 02, which lead to three physical phases. The PMNS matrix in Eq. (IBaj) 
can, equivalently, also be directly written as 



PMNS 



f/-diag(e 1 ^ 1 ,e 102 ,l) 



(8b) 



where U is a CKM-like matrix that is on the form as in Eq. (15bl) and the phases 0i and 02 
are Majorana phases. The CKM-like matrix U in Eq. (18b|) is described by the solar angle 
6*i2, the reactor angle 6*13, the atmospheric angle 6*23, and one Dirac CP-phase 5, which we 



identify in the standard parameterization of Eq. (15bl) as 9 



9ij and S 



5. The PMNS 



matrix has thus 3 mixing angles and 3 phases and contains therefore 6 physical parameters. 
Writing the matrix elements of C/pmns as = (£7pMNs)ij> we rea d off in the standard 
parameterization of Eq. f!5bj) . the 6 leptonic mixing parameters as follows 



sin 6>i3 
tan 6*12 



U 13 \ 

U 12 



tan 9 



23 



Uxx 
U 23 



h2 = - if Uxx 



33 



723 



~ if ^33 = 0, 



arg 



^12^13^22^23 



-Sl2"Sl3S23Ci2C 13 C23 C12C23 / 



+ a J*™«\ for 0,|, 



(9a) 
(9b) 

(9c) 
(9d) 



else, 



arg(e- l5 t/* 3 f/ 11 ), 

arg(C/ 31 C/ 3 *3), 



arg(e- iS U^U 



13^12; 



arg^C/*^), 

arg(C/ 22 C/ 2 * 3 ), 





for 9 {j ^ f and 13 ^ 0, 
for 9ij 7^ | and only #13 = allowed 
else, 

for 9 tj I and 12 ,0 13 ^ 0, 
for 0i3,e 2 3 = O and <9 12 ^ 0, f 5 
for or = but + f , and 9 23 ^ 
else. 



(9e) 



(9f) 



' 2' 



Note that we set the phases to if they are undefined. Since we can easily identify these 
cases if needed, there is no bias or constraint introduced by this choice. For example, for 
6*i3 = 0, 5 is undefined. In this case, 5 does not affect the mixing matrix, i.e., an arbitrary 
choice of does not change the physics. However, for the phase prediction of S, the irrelevant 
cases can be easily eliminated by the identification of the corresponding #13 = 0. Note that 
the above relations are valid for a general unitary matrix £/ un i ta ry as in Eq. (I5a| as well. 
In addition, if £/pmns is already a CKM-like matrix, i.e., all unphysical phases are already 
removed, we simply have 5 = arg([/* 3 ). 



5 



2.2 Generating the PMNS Matrices and Textures 



Let us now introduce our method for a systematic scan of the textures that exhibit QLC 
Instead of starting out with various forms for the lepton mass matrix textures and calculate 
from these the lepton masses and mixing angles, we follow in our procedure a reverse 
process: we reconstruct the mass matrix textures in flavor basis from the masses and mixing 
parameters of charged leptons and neutrinos. 

Although neutrino oscillation experiments will in future allow to pin down the leptonic 
mixing angles with increasing precision, the PMNS matrix does not uniquely reveal the 
individual mixing parameters in Ug and U v . Existing models and studies mostly suppose that 
the observed large solar and atmospheric mixing angles arise mainly in the neutrino sector. 
It is, however, important to emphasize that maximal atmospheric mixing can also equally 
well arise from the charged lepton sector as proposed, e.g., in "lopsided" GUT models [39] 
(for realizations of bilarge neutrino mixing in lopsided models see also, e.g., Refs. [40-42]). 
Moreover, the QLC relations in Eq. (j4]) suggest, in particular, that both maximal (~ vr/4) 
and small (~ 0c,0 c b) mixing angles might be expected in the charged lepton and in the 
neutrino sector. Motivated by QLC, we will thus assume in our approach that the mixing 
angles 0f • in Ug and Q\a in U v can a priori take any of the values in the sequence 7r/4, e, e 2 , . . . , 
where e ~ 9q ~ 0.2, and then compare the resulting PMNS mixing angles 6ij with current 
data. The choice of the angles in this sequence is a simple and straightforward interpretation 
of the QLC relations in Eq. in the sense that the solar angle 6* 12 — 33° can only arise as 
a result of taking the product of Ug and U v in ?7pmns = U\U V . The assumption that any of 
the angles 0f ■ and 9^ can assume any of the values 7r/4, e, e 2 , . . . , generalizes the definition 
of QLC to what we call an "extended QLC". 

Let us next specify the values of the mass parameters in our approach. Motivated by the 
mass ratios of quarks and leptons in Eqs. (T5|) and (j3J), we will approximate the diagonalized 
charged lepton mass matrix M diag in Eq. flTaj) by 



and choose for the diagonalized neutrino mass matrix M^ iag in Eq. ( 17bl) any of the cases 



= m 3 diag(e 2 , e, 1), = m 2 diag(l, 1, e), = m 3 diag(l, 1, 1), (10b) 



where the first, second, and third equation corresponds to a normal, inverted, and degenerate 
neutrino mass spectrum, respectively. We thus see that the quantity e ~ 9q will appear in 
our approach as a small expansion parameter, which is in charge of all the hierarchies and 
small effects that are exhibited by the mass spectra and mixing angles of the fermions. 

Our systematic search for textures satisfying QLC hence consists in the following simple 
three-step procedure: 

First step - In the first step of our procedure, we consider the PMNS matrix in the 
parameterization £/pmns = U\U V = U\DU V K of Eq. ( J5al) and assume here that the entries 
sin(0|) and sin(6^) in the CKM-like matrices Ug and U v can take all the values 



Ml 



diag 



m T diag(e , e 




(10a) 



4' 4 e {^p e > e2 >°}> 



(11) 
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where we have defined sf 3 - = sin(0f-) and = sin (6^). For definiteness, we choose in 
Eq. (fTTj) e = 0.2 as our standard value. Moreover, we will first restrict ourselves to the case 
of real mixing matrices Ue and U v , and assume in Eq. (IBa| that the phases S e , 8 V , tpi, tp2, 4>i, 
and 02, are only or n. In total, this gives 



possible distinct pairs of mixing matrices {Ue,U v } (cf. Eq. (IHal) ). Later, in Sees. 14.21 and 
I4.3[ we will vary e in the range 0.15 < e < 0.25 and the phases S e and S u in the whole range 
8 l , S u G [0, 271"]. In the following, we will call a "model" the set of twelve mixing parameters 

{^12; ^13 > ^23; ^12 > ^13 j ^23; " j ^ V i <Pi> <f2j 4>\i 02} that describes a pair of charged lepton and 
neutrino mixing matrices {Ue,U u }. 

Second step - In the second step of our method, we obtain for each of the 262 144 matrix 
pairs the corresponding PMNS matrix {Jpmns = U\U V , read off the PMNS mixing parameters 
012; 013, 023, 5, 0i, and 02 using Eq. (jUj), and extract those pairs {Ue,U u }, which are in 
agreement with current neutrino oscillation data. In order to have an automatic selection 
criterion for our candidate pairs {Ue, U u }, we create our "best sample" by defining a selector 



This selector corresponds to a Gaussian x 2 with the current best-fit values and the given 
relative la errors cr 12 ~ 9% (for sin 2 #12) and 023 — 16% (for sin 2 #23) from Tabled], where 
we assume a Gaussian distribution in sin 2 #12 and sin 2 #23 as an approximation. For the best 
sample, we choose all models that satisfy the selection criterion 



which corresponds to a A% 2 for the 3a confidence level with two degrees of freedom, and the 
hard cut on #13 represents the current 3a bound according to Table [U We do not include 
the numerical value of sin 2 #13 in the selector in Eq. fTTSl . because its current best-fit value 
would introduce a bias towards small #13. In some cases, we will use different values for 
the selection, such as to test the experimental pressure on the model space. For example, 
on a time scale of about ten years from now, we choose the values from the last column in 
Table [T] This selection process is a conservative guess/crude estimator for the sample of 
models which can be still accommodated with data within the 3a confidence level. Note 
that we have, so far, neither used the mass squared differences (and mass hierarchy), nor 
extracted/predicted them. At this stage, in steps 1 and 2, our procedure is independent 
from the mass spectra of charged leptons and neutrinos. They become, however, important 
for constructing the mass matrix textures in the next step. 

Third step - The third step of our approach consists of considering all the pairs {Ue, U v } 
that have been selected in the previous step 2 and rotate for each such pair the diagonal 
matrices Mf ag and ag given in Eqs. (flOl) back to flavor basis according to 



4 3 x 2 x 4 3 x 2 5 = 262 144 



(12) 




(13) 



S < 11.83 and sin 



2 013 < <Ti3 ^ 0.04, 



(14) 



M e = UeM° 



diag 



M, Maj = U v Mf^U T v , M\ 



Dirac 



U V M\ 



diag 



(15) 



V 



V 
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where we have chosen, for simplicity, U' e = U' u = I3. The Dirac mass matrices Me and 
j^Dirac can be viewed as representatives of a class of mass matrices that is obtained 
by introducing arbitrary rotation matrices U' e and U' u acting on the right-handed leptons. 
Moreover, note that Eq. (1151) actually describes six different cases, depending on whether one 
chooses from Eq. (II Obi) a normal hierarchical, inverted hierarchical, or a quasi-degenerate 
neutrino mass spectrum. In the following, we will, as already mentioned in the introduction, 
denote the mass matrices Me, M^ 1 ^, and M® nac , in Eq. (fT5j) expanded in terms of e as "mass 
matrix textures" or simply "textures" . 

Let us briefly summarize again our three-step procedure for a systematic scan of real PMNS 
matrices and textures with QLC: 

1. Generate 226 144 PMNS matrices by inserting all distinct combinations of sf •, s^-, 5 , 

8 W , (pi, (p 2 , 0i, and 02, into Eq. flHal) . Here, sfj and s?- are taken from {-^,e, e 2 ,0}, 
whereas the phases can assume values 7r and 0. 

2. Read off the mixing parameters #12, #13, #23, 5, <pi, and 4>2, using Eq. (jUJ) and select 
the pairs {Ue, U u } where C/pmns = U\U V satisfies the selection criterion in Eq. (fT4l) . 

3. For each selected pair {Ue, U u } rotate the mass matrices in Eq. (fTUl) back to flavor 
basis to obtain the textures M e , M^ aj , and Mf ac (cf. Eq. (fT5]) j. 

We wish to emphasize that the above procedure leaves complete freedom to let the observed 
leptonic mixing angles be generated by contributions from both the charged lepton and/or 
neutrino sector in the most general way. Another advantage of our method, that we will 
exploit later, is that it allows to scan very quickly the parameter space of complex phases by 
varying the parameters directly in the PMNS matrix instead of varying the highly redundant 
phases in the mass matrices. Our procedure to construct mass matrix textures in this way 
is extremely simple and efficient. In particular, to derive the textures, there is no need to 
perform any diagonalization in the whole process. 

3 Scanning the Parameter Space for Individual Models 

In this section, we apply our procedure outlined in Sec. 12.21 for extracting the mass matrix 
textures from the best sample obtained with the selector in Eq. ({TBI) , and make it explicit 
for several examples. Thereby, as explained above in step 3 of our procedure, we start with 
the diagonal forms of the charged lepton and normal hierarchical neutrino mass matrices in 
Eqs. ( II 01) and obtain the mass matrices by rotating back to flavor basis following Eq. (I15p . 
where we set for convenience U[ and U' v equal to the unit matrix. In the following, we make 
first a subjective choice and focus in Sec. l3.1l on several textures which are of special interest 
like (i) the combination of a CKM-like matrix for Ue with a bimaximal mixing matrix U v , 
(ii) a texture with very small #13, and (iii) the sets which represent the best approximation 
to the current best-fit values, i.e., those for which the selector S in Eq. ( TT3"j) takes the lowest 
value. Later, in Sec. 13.21 we present a complete list of the textures that exhibit a small #13, 
as well as survive increased experimental pressure from the other parameters. 
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3.1 Several Examples of Models Compatible with Current Bounds 



Here, we present a subjective selection of models compatible with current bounds. In 
addition, we describe how we identify the leading order terms in the textures. 



Matrices with CKM-like Plus Bimaximal Mixing 



Let us first consider the case where Ui is CKM-like and U u is on a bimaximal mixing form. 
In the parameterization of Eq. fl8a|) . this can be realized in our procedure by taking for the 
charged lepton and neutrino mixing parameters 



V <f 



0l2, S 13> S 23> <Pl, $2, 01, 02) 

which lead to the PMNS mixing angles 



(e,0,e 2 ,vr) 

11 

Vf' 6 ' '°' '°' 7r ' 7r ' ) ' 
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The corresponding charged lepton, neutrino, and PMNS mixing matrices are 
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0.8 e 2 -0.5 + 1.2 e + 0.8 e 2 
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(17a) 



(17b) 



The CKM matrix and the bimaximal mixing matrix are, on the other hand, given by 



Vc 
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and U] 
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2 
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L 1_ 

2 72 

L J_ 

2 v/2 



(17c) 



where we have used for Vckm the Wolfenstein parameterization to second order in A = 

0.22. The PMNS mixing angles in t/pMNS — VcKM^bimax are (#12, $13, $23)cKM+bimax — 

(36°, 9°, 45°). Comparison of Eqs. fll7al) and fU7cl) yields that Ue reproduces exactly the 
CKM matrix in the Wolfenstein parameterization, whereas the entries of U v differ from 
those in L^max by terms of the orders e or e 2 . These terms lead to the deviations between 
the numerical values of (6 I 12 , #13, #23)cKM+bimax and (^12,^13,^23) m Eq. (j 16 eft . Note that 
choosing the phases of U v such that all signs of the entries in f/bimax are reproduced but 
without changing Ue, changes the PMNS mixing angles in Eq. (1 16c[) to (#12, #13, #23) — * 
(33.8°, 10.0°, 43.4°). 

Now, we rotate the mass matrices M^ iag = m T diag(e 4 , e 2 , 1) and M^ iag = m 3 diag(e 2 , e, 1) 
back to flavor basis following Eq. (fl5]) . where Ui and U v are parameterized in terms of 
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powers of e as in Eq. (1 1 7aj) . In flavor basis, the resulting mass matrix textures of charged 
leptons and neutrinos are then given by 
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2^ 
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where we have, for simplicity, assumed the special case of a normal hierarchical neutrino 
mass spectrum. Starting, instead with M^ iag = diag(l, 1, e) or M^ iag = m 3 diag(l, 1,1), 
we obtain the corresponding textures of the neutrinos for the case of inverse hierarchical or 
degenerate neutrino masses. In Eqs. (JTB1) . "— symbolizes, up to an overall mass scale, the 
identification of the leading order terms in the expansion in e that contribute to the mass 
matrix elements in Me, M^ 1 ^, and M° irac . This identification is possible for all the mass 
matrices which are currently valid at the 3a level, since it turns out that the contributions 
to a mass matrix element belonging to different orders in e can always be clearly separated 
from each other. In other words, for the given expansions of the mass matrix elements in 
powers of e, the zeroth order term in e (if non-zero), is always larger than the higher order 
terms, and the first order term in e (if non-zero) is always larger than the second order term. 
In the thus identified leading order contributions to the mass matrices, we can then further 
approximate and set the order unity coefficients equal to one. This leads us finally to a rough 
texture with matrix elements 1, e, e 2 , and 0. These entries have thus to be understood as 
order of magnitude entries. Note that such textures become 't Hooft natural [43] when they 
arise from a spontaneously broken flavor symmetry [44] (see also Ref. [45]). An origin of the 
small number e is provided in terms of an anomalous U(l) symmetry [46], which has been 
employed for generating textures in various models [47] (for an anomaly-free approach see, 
e.g., Ref. [48]). 

In the rest of the text, we will always proceed exactly the same way and determine the 
textures for charged leptons and neutrinos as it was done above in arriving at Eqs. (1181) . 



Matrices with Very Small 6 1S — Near-CKM Plus Bimaximal Mixing 

We find a class of 64 models with very small #13 ~ 0.2°. All these models have in common 
that 9 12 is slightly off the current best-fit value, which means that they will be excluded 
in ten years by our assumptions. Note, however, that experimentally 9 13 and # 12 are two 
different degrees of freedom, which means that there could well be much stronger pressure 
coming from #13 than from 9u- In this case, this class of models may be surviving very 
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long. For one typical representative combining near- CKM- type mixing in Ui with with 
near-bimaximal mixing in U u , the charged lepton and neutrino mixing parameters are 
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(19a) 
(19b) 



which lead to the PMNS mixing angles 

(012,013,023) = (28.7°, 0.2°, 46.1°). 
The charged lepton, neutrino, and PMNS mixing matrices are 
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(21a) 



(21b) 



4 2V2 / 

Here, C/^ is on a bimaximal mixing form but U v is not exactly the CKM-matrix (mainly 
due to the entry ~ e in the 1-3 and 3-1 element). For the textures of charged leptons and 
neutrinos we thus find 
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which is similar to the example given before. 





(21d) 



(21e) 



Best-Fit Matrices 

Let us now consider the two best examples of models which provide perfect fits to current 
data. These two examples minimize the selector in Eq. (|T3|) . which takes for these models 
the smallest possible value S = 0.12. We denote these textures as "best-fit matrices" since 
they provide a "perfect fit" to current data. Both examples exhibit similar PMNS mixing 
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angles with a 6*13 that is in the reach of next-generation neutrino oscillation experiments. 
For example, 6*13 is large enough to be discovered by Double Chooz, which means that for 
large #13, these representatives would be the perfect candidates. 

Example (a) - Our first example for the best-fit matrices has the charged lepton and 
neutrino mixing parameters 



1*125 *13> ^23? ) 
( S 12> S 13> S 23> £l» ^2, 01, 02) 

which lead to the PMNS mixing angles 
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(22b) 



(012,013,023) = (33.4°, 7.5°, 43.5°). 



(23) 



The charged lepton, neutrino, and PMNS mixing matrices are 
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(24a) 
(24b) 



Note that, in this case, JJg is on a bimaximal mixing form but U u is neither the CKM-matrix 
nor CKM-like. For the textures of charged leptons and neutrinos we find 
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(24e) 



We thus see that Mi is on a typical lopsided form with a strictly hierarchical structure 
within each row. On the other hand, M^ Ia -' and M i I p irac have nearly degenerate large entries 
in the first and third row/column with small entries in the second row/column. 
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Example (b) - Our second example for the best-fit matrices has the charged lepton and 
neutrino mixing parameters 



'12, °13' °23 
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which lead to the PMNS mixing angles 

(012,013,023) = (33.4°, 7.5°, 43.5°). 
Thus, the charged lepton, neutrino, and PMNS mixing matrices are 
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As a consequence, we find for the textures of charged leptons and neutrinos 
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(25a) 
(25b) 

(26) 



(27a) 



(27b) 



(28a) 



(28b) 



(28c) 



In this case, Mg is no longer of the usual lopsided form but has only entries of similar (small) 
orders in the first and second row. The neutrino mass matrices, however, have entries in 
the last column that differ only by order unity factors, which leads, after neglecting these 
factors, for M^ aj to a "democratic" or "anarchic" texture. 
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3.2 Textures with Small 6 13 Surviving Increased Experimental Pressure 

Let us now concentrate on models with a small #13 that would survive a projected increased 
experimental pressure in ten years, i.e., pressure from a stronger # 13 bound as well as from 
the other parameters as estimated in the last column of Table [TJ The reason for considering 
textures with small #13 is two-fold. First, from a theoretical point of view, an understanding 
of textures with small #13 seems to require interesting model-building assumptions, e.g., 
discrete non-Abelian flavor symmetries [49]. Second, from an experimental perspective, a 
general survey shows that most of the existing GUT model-building approaches predict a 
comparatively large 6*13, which is close to the current upper bound [50]. These predictions 
from model-building can soon be tested in next-generation experiments such as Double 
Chooz, i.e., the measurement of #13 exerts a strong experimental pressure on currently 
allowed models. 

In Table [21 we have summarized all possible real textures (only leading order entries) that 
would survive a projection of the experimental bounds on the leptonic mixing angles in ten 
years. In order to extract these, we have used the extrapolated errors from the last column 
in Table [T] in Eq. (TT51) and Eq. <HM by assuming that the current best-fit values hold. Of 
course, we could as well extract the textures for any other set of best-fit values with the 
same procedure, but this would exceed the scope of this paper (see also Ref. [60]). We have 
listed the textures Me, M^ a \ and MjP irac , for normal hierarchical, inverted hierarchical, and 
degenerate neutrino masses. The textures in Table [2] have been selected from the 264 144 
initial models following the procedure in Sec. 12.21 which assumes that the charged lepton and 
neutrino mixing matrices Ue and U v are real. The textures in Table [2] therefore represent 
all the models which survive the projected experimental bounds in ten years for the CP 
conserving case. The number of such models is less than ~ 0.01% of the initial sample. 

In the notation and parameterization of Sec. 12.21 Table [2] shows the individual mixing 
parameters (s[ 2 , s{ 3 , s 23 ) , (s\ 2 , £23)) an d (S e , 5 V , (pi, <p 2 ), of Ue and U v , as well as the 
PMNS mixing angles (#12, #13, #23)- Notice that we have dropped the phases (pi and 02, 
since they can be absorbed into the definition of the Majorana phases <p\ and <p 2 . In the 
CP conserving case considered here, however, the Majorana phases are either or tt, and 
thus play no role for the structure of the textures (or the PMNS mixing angles). Moreover, 
Table [2] contains for degenerate neutrino masses only the Dirac neutrino mass matrix M^ irac , 
since in the CP conserving case M^ 1 ^ would only be proportional to the unit matrix. 

For the textures #6, #11, and #12, we have given two sets of PMNS mixing angles, which 
correspond to different choices of the phases and n in Ue and U u , resulting in a variation 
of the mixing parameters. This is not the case for all other textures in the list. Moreover, 
even though different phase combinations change the PMNS mixing angles, the textures of 
the mass matrices are, for normal hierarchical neutrino masses, invariant under such permu- 
tations of the phases. For an inverted hierarchical or degenerate neutrino mass spectrum, 
however, a different choice of phase combinations can manifest itself in a slightly different 
mass matrix texture, even in the CP conserving case. 



14 



# 




Normal ] 

M Maj 


iierarchy 

M Dirac 


Inverted ] 

M Maj 


iierarchy 

M Dirac 


Degenerate 

M Dirac 


( s 12' s 13' s 23) 
( s 12i s 13> s 23) 


(#12, #13> #23) 


1 




/ e \ 

e 2 
\ 1 / 


fee e \ 

ell 
\f 1 1/ 


/ e 2 e e \ 

e 2 e 1 
V e 2 e 1 / 




(iii) 




( 


( 1 1 e 2 \ 

1 1 e 
V 1 1 e j 


I 




I I e\ 
111 

V 1 1 i/ 


(e 2 , e, 0) 

(£, TT, 0, £ + 7r) 


(35.2°, 4.9°, 43.8°) 


2 


1 


( e \ 

e 2 e 2 

V 1 / 


fee e \ 

ell 
\f 1 1/ 


/ e 2 e e \ 

e 2 e 1 
\ e 2 e 1 / 




Cii) 




I 


( 1 1 e 2 \ 

1 1 e 
V 1 1 e , 


I 




/I 1 e\ 
111 

V 1 1 i/ 


(tt, tt, 0, 0) 


(35.5°, 4.5°, 41.6°) 


3 


1 


/ e 2 \ 
e 2 

V 1 / 


fee e \ 

ell 
\f 1 1/ 


/ e 2 e e \ 

e 2 e 1 
V e 2 e 1 / 








( 


( 1 1 e 2 \ 

1 1 e 
V 1 1 e , 


I 




/I 1 e\ 
111 

V 1 1 i/ 


(e, e 2 , 0) 
(£, 0, 0,£ + tt) 


(35.2°, 4.8°, 46.6°) 


4 


1 


( e 2 \ 

e 2 e 2 

V 1 / 


let e \ 

ell 
\f 1 1/ 


/ e 2 e e \ 

e 2 e 1 
V e 2 e 1 / 




Cii) 




I 


M 1 e 2 \ 

1 1 e 
V 1 1 e , 


I 




/I 1 e\ 
111 

\i 1 i/ 


(e, e 2 , e 2 ) 

(72' e ' 75^ 
(0, 0, 0, tt) 


(35.5°, 4.5°, 48.9°) 


5 




/ e \ 

e 2 
\ 1 / 


fee e \ 

ell 
\ e 1 1 / 


/ e 2 e e 2 \ 

e 2 e 1 
V e 2 e 1 / 


( 


1 e 2 e 2 
e 2 1 1 
k e 2 1 1 


) 




(!!:) 




( 


/ 1 1 e 2 \ 

111 
,111/ 


(e 2 , e, 0) 
(±- e 2 -J-l 

(£, TT, 0,£ + 7r) 


(35.2°, 3.7°, 45.1°) 


6 


1 


/ e \ 

e 2 e 2 

V 1 / 


fee e \ 

ell 
\ e 1 1 / 


/ e 2 e e 2 \ 

e 2 e 1 
V e 2 e 1 / 


( 


' 1 e 2 e 2 
e 2 1 1 
e 2 1 1 


) 








( 


< 1 1 e 2 \ 
111 

v i i i ; 


(c 2 , e, e 2 ) 

(0, TT, 0, TT) 

(tt, tt, 0, 0) 


(35.0°, 3.8°, 47.7°) 
(35.5°, 4.6°, 43.1°) 



# 




Normal ] 

M Maj 


iierarchy 

M Dirac 


Inverted ] 

M Maj 


iierarchy 

M Dirac 


Degenerate 

M Dirac 


( s 12' s 13' s 23) 
( s 12i s 13> s 23) 

(6^5^,^,02) 


(#12) #13> #23) 


7 


\ 


'00 e \ 

e 2 
v 1 / 




(Mi) 






/ 2 2 \ 

(MM 




I 


/ 1 2 2 \ 

/ 1 e^ e 2 > 

e 2 1 1 
v e 2 1 1 y 


1 




(MM 






/ 1 1 2 \ 

/ 1 1 e z \ 

iii 
V i i i J 


(e, e 2 , 0) 
e 2 -M 

K, o, o,e + ^) 


(35.2°, 3.7°, 45.1°) 


8 


I 


f e 2 \ 

e 2 e 2 
v 1 / 




(Mi) 






(MM 




I 


/ 1 e 2 e 2 N 

e 2 1 1 
V e 2 1 1 J 


I 




(MM 






/ 1 1 e 2 \ 

V i i i / 


(e, 6 2 , e 2 ) 

(tt, 0, 0, 0) 
(0, 0, 0, tt) 


(35.0°, 3.9°, 42.8°) 
(35.5°, 4.7°, 47.3°) 


9 




/ e \ 

e 2 1 
V e 2 1 / 


1 


(e e e \ 
e e e 2 

V e e 2 i y 


I 


( 


' e 2 e e N 

e 2 e e 2 
v , 2 1 y 


) 








I 


( 1 1 e 2 \ 

1 1 
y e e e ; 


1 




/ 1 1 e \ 

lie 2 
V e e 1 ; 


( e2 > e ' 71 ) 

e < e2 ) 
(0, 0, 0, 0) 


(35.5°, 4.5°, 41.6°) 


10 


( 


/ e 2 \ 

e 2 1 
V e 2 1 / 


1 


/e e e \ 
e e e 2 

V e e 2 i y 


1 


( 


' e 2 e e ^ 

e 2 e e 2 
V , 2 1 , 


) 




Ok) 




I 


n i £ M 

1 1 

y e e e J 


1 




/ 1 1 e \ 

lie 2 
\ e e 1 / 


(0, 0, 0, tt) 


(35.5°, 4.5°, 48.9°) 


11a 




/ e \ 
e 2 1 

\ e ' 1 / 


( 


' e e e 2 " 

e e e 2 
v e 2 e 2 1 


) 




(:);<) 




1 


(l e 2 \ 

1 e 2 
V e 2 e 2 e , 


) 


I 


( 1 1 \ 
1 1 

v e 2 ey 


1 


( 


' 1 1 e 2 \ 

lie 2 
v e 2 1 / 


e ' 75) 

(0, TT, 0, TT) 


(35.5°, 4.7°, 47.3°) 


lib 




/ e \ 

e 2 1 
V e 2 1 / 


( 


' e e e 2 " 

e e e 2 
v e 2 e 2 1 , 


) 




(::::)) 




I 


/ 1 e 2 \ 

1 e 2 
V e 2 e 2 e ) 


1 


I 


f 1 1 \ 
1 1 

V e 2 e J 


1 


( 


' 1 1 e 2 \ 

lie 2 
v e 2 1 / 


75) 

f^- e 2 e 2 ) 
(0, 0, 0, 0) 


(35.0°, 3.9°, 42.8°) 



# 




Normal ] 

M Maj 


iierarchy 

M Dirac 


Inverted ] 

M Maj 


iierarchy 

M Dirac 


Degenerate 

M Dirac 


( s 12' s 13' s 23) 
( s 12i s 13> s 23) 

(6^5^,^,02) 


(^12) #13> ^23) 


12a 


( 


( e 2 \ 

e 2 1 
V e 2 1 / 


/ e e e 2 \ 

e e e 2 
V e 2 e 2 1 / 




/ e 2 e e 2 \ 
e 2 e e 2 ) 

V 1 / 


( 


' 1 e 2 ' 

1 e 2 
v e 2 e 2 


) 




K l e 2 lj 




( 


( 1 1 e 2 N 

1 1 e 2 
yO e 2 1 , 


1 


<?> 7f ) 

(0, TT, 0, 0) 


(35.5°, 4.6°, 43.1°) 


12b 


I 


^00 e 2 \ 
e 2 1 

v e2 1 / 


/ e e e 2 \ 

e e e 2 
V e 2 e 2 1 / 




^ e 2 e e 2 \ 
e 2 e e 2 \ 

V 1 / 


( 


' 1 e 2 N 
1 e 2 


) 




J 2 o % 




1 


(I 1 e 2 N 

1 1 e 2 
v , 2 1 , 


1 


e2 > 75) 
f^- e 2 e 2 ) 

(0, 0, 0, tt) 


(35.0°, 3.8°, 47.7°) 


13 




/ e \ 

e 2 1 
V e 2 1 / 


fee e \ 

e e e 2 
V e e 2 1 / 




/ e 2 e e \ 

e 2 e 
V e 2 1 J 




(iSi) 






y e e e J 






(Iii) 




e ' 75) 

(o, 0, 


(35.2°, 4.9°, 43.8°) 


14 


( 


/ e 2 \ 

e 2 1 
V e 2 1 / 


fee e \ 

e e e 2 
V e e 2 1 / 




( e 2 e e\ 

e 2 e 
V e 2 1 / 




(iii) 






(i i •) 






(Iii) 




e2 > 75) 
(73> e ' °) 

(0, 0,^ + Tr) 


(35.2°, 4.8°, 46.6°) 


15 




/ e \ 

e 2 1 
V e 2 1 J 


/ e e e 2 \ 

e e 
V e 2 1 / 




( e 2 e e 2 \ 
e 2 e 

V 1 / 


1 


f 1 e 2 \ 
1 

v e 2 o e y 




( 


' 1 1 , * 
1 1 

v e 2 e 2 e y 


) 


( 


' 1 1 e 2 ' 

1 1 
v £ 2 £ 2 1 , 


) 


e ' 71) 
e2 ' °) 

(o, 0, 


(35.2°, 4.2°, 45.4°) 


16 


( 


{00 e 2 \ 

e 2 1 
yO e 2 1 J 


( e e e 2 \ 

e e 
V e 2 1 / 




^ e 2 e e 2 \ 
e 2 e 

V 1 / 


1 


M e 2 \ 

1 
V e 2 e ) 


i 


( 


' 1 1 N 
1 1 

,e 2 e 2 e, 


) 


( 


' 1 1 e 2 ' 

1 1 
v £ 2 £ 2 1 , 


) 


e2 > 75) 
e2 ' °) 

(0, 0, e + 


(35.2°, 3.7°, 45.1°) 



# 


M, 


Normal ) 

M Maj 


iierarchy 

?\#Dirac 


Inverted ] 

M Maj 


iierarchy 

i\/T Dirac 

V 


Degenerate 

;\/rDirac 

V 


( S 12i S 13> S 23) 
( S 12) s 13' S 23) 


(#12, #13i ^23) 


17 




[Hi) 




( 


' 1 e 2 1 x 
e 2 e e 2 

v 1 e " 1 


) 




( £2 t2 M 

e 






(iii) 




1 


e 1 

v 1 « «y 






(iii) 




^75' e ' 75' 1 
( e ' o) 


(35.2°, 3.8°, 50.8°) 


18 








( 


' 1 e 2 1 x 
e 2 e e 2 

v 1 e " 1 


) 


( 


' e 2 e 2 1 > 

e e 2 
v £ 2 e 2 1 , 


) 


( 


1 e 2 1 
e 2 1 e 2 
1 e 2 1 


) 


l 


e 1 




( 


' 1 e 1 N 

e 1 e 2 
v 1 e 1 , 


) 


f— e 4=1 

(7T, 7T, 7T, 0) 


(33.6°, 3.1°, 52.2°) 


19 




(ifi) 






(ill) 






(si'i) 






(Iil) 






(iii) 






(iii) 




e ' °) 

fe -J- 4=1 


(35.2°, 3.8°, 50.8°) 


20 


( 


' ( s e ^ 

« 2 e 2 
,001, 


) 




(iii) 












(iii) 






(iii) 






(iii) 




«, ^) 

( e > 73' 73^ 

(7T, 7T, 7T, 0) 


(33.6°, 3.1°, 52.2°) 



Table 2: Mass matrix textures surviving the projected increased experimental pressure in ten years (see Table [T]). For the mass matrices of 
charged leptons and Dirac neutrinos (for normal hierarchical, inverted hierarchical, and degenerate neutrinos, respectively), we use the basis in 
which the mixing angles of the right-handed fields are zero. Note that in the case of degenerate neutrinos, the Majorana mass matrix textures 
are simply the unit matrix and they are therefore neglected. The corresponding mixing angles and phases of Ue, U v (with £ € {0, 7r}) as well as 
the mixing angles of E/pmns are a l so shown. Note that in some cases, a phase variation in Ui and J7„, without changing their mixing angles, can 
lead to different forms of E/pmns- I n these cases, we relate the causal phases of Ug and U u to their corresponding PMNS mixing angles. 



The textures in Table [2] can be roughly divided into four classes: 

1. "Semi-anarchic" textures (#1-8): Q V X2 and #23 are maximal, all other mixing angles 
are small. Here, Mg has a strongly hierarchical structure, whereas M^ Ia -' and M° irac 
have several leading order entries. 

2. Lopsided textures (#9-16): #23 an d #12 are maximal, all other mixing angles are small. 
Mg is on a lopsided (highly asymmetric) form. 

3. "Diamond" textures (# 17-18): 8[ 2 , #23> an d #13, ar e maximal, all other mixing angles 
are small. Mg is on a lopsided form. In the neutrino mass matrix, leading order entries 
are situated in the corners. 

4. Anarchic textures (#19-20): 6[ 2 , 6^, and 8 23 , are maximal, all other mixing angles 
are small. Mg has a strongly hierarchical structure, whereas M„* and M^ lTac have 
many leading order entries. 

We have called the matrices #17 and #18 "diamond" textures, because of the diamond-like 
structure of the small entries in . Most of the textures in Table [2] fall either into the 
first or the second class: 8 textures are semi-anarchic and 8 are lopsided. The remaining 
four textures are equally distributed between the diamond and the anarchic forms. For 
the semi-anarchic and lopsided textures, the charged lepton and neutrino mixing matrices 
exhibit each a single maximal mixing angle. The diamond and anarchic textures, on the 
other hand, have in Ug (diamond) or U v (anarchic) an additional maximal mixing angle such 
that the total number of maximal mixing angles in the charged lepton and neutrino sectors 
is three. 

The diamond textures (#17 and #18) are interesting, since they can be viewed as represent- 
ing a generalization of the lopsided models with two maximal mixing angles in the charged 
lepton sector. They also have (like the anarchic examples) the unusual feature that 6i 3 is 
maximal. Let us therefore briefly consider these mass matrices in some more detail. For 
the diamond textures, the charged lepton, neutrino, and PMNS mixing matrices are 
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We thus obtain for the expansion of the diamond textures in powers of e the forms 
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Keeping in the expansion in e only the leading order coefficients and approximating these by 
one, we arrive from Eqs. (!3~T!) at the textures #17 and #18 in Table El The best-fit matrices 
of Sec. 13 .11 that provide a perfect fit to current data are not contained in Table [2J since they 
exhibit a reactor angle #13 = 7.5° that would be ruled out by the projected confidence levels 
in 10 years. The diamond textures #17 and #18 can, however, be considered as a crude 
approximation of the best-fit matrices. 

It is interesting to compare the results in Table [2] also with tri-bimaximal mixing [51], 
which has the PMNS angles (0i 2 , #13, 2 3)tribi — (35°, 0°, 45°). In our procedure, neither Ui 
or U v are on a tri-bimaximal form, since our parameter space does not contain s^ r 2 lbl = 
1/V3. Nevertheless, we can associate some textures in Table [2] with perturbations to tri- 
bimaximal mixing. For instance, the textures #5, #7, #15, and #16, could be viewed 
as describing (small) deviations from tri-bimaximal mixing (deviations from tri-bimaximal 
mixing have been considered previously, e.g., in Ref. [52] and references therein). Different 
from the literature, however, none of the examples for tri-bimaximal-like mixing in Table [2] 
is associated with a charged lepton mixing matrix matrix U( that is close to unity. The 
tri-bimaximal-like mixing emerges in these cases, instead, from large mixings in both the 
charged lepton and the neutrino sector. 

Up to corrections of the order ~ e 3 , all the data models (in the sense of Sec. 12. 2j) for the 
textures #1-16 in Table [2] exhibit values for $12 that follow the relation 

012 + 4= + 4= = 7" ( 32 ) 

This corresponds to the usual QLC relationship in Eq. (j4j). The models for the diamond 
and anarchic textures #17-20 satisfy, instead, the new sum rule 

3 

012 + 7= e = arctan(2 - y/2). (33) 

5 + 2y2 

This new relation is an outcome of our extended QLC approach. The models for the textures 
#1-16 exhibit 

013 = 0(e) and 23 = j + 0(e 2 ), (34) 
whereas the models for the diamond and anarchic textures #17-20 reveal 

013 = arcsin(i(2 - y/2)) . 1 e , (35a) 

4 V5 + 2y/2 

23 = arctan(l + —=) + —(2 - Uy/2) e . (35b) 
V 2 17 
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In Eq. (I34p . #13 is proportional to e, but it is nevertheless small due to a small coefficient 
multiplying e. In contrast to this, in Eq. (l35aD . a small #13 arises from ~ e corrections to a 
small 0(1) ~ 8° term. A close to maximal #23 results in Eq. (135b[) from a ~ e correction 
to arctan (l + l/y/2) ~ 60°. We would like to emphasize again that the new sum rules in 
Eqs. (133]) and fl35l) are, together with the known relations of Eqs. (152]) and (134"]) . consequences 
of the extended QLC approach. 

4 Distribution of Observables: Mixing Angles 

In this section, we investigate the ensemble of models found in our best sample. This includes 
all models as defined by our selector (cf. Eq. (TIB"]) ) for the currently allowed parameter 
ranges. We will be mainly interested in distributions of observables from this ensemble of 
models on a statistical basis rather than individual models. Note that there will be no 
predictions for the mass hierarchy, as well as there will be no dependence of the predictions 
on the hierarchy. The reason for this is the independence of the selector and the model 
building process on the mass hierarchy (only the neutrino mass textures depend on the 
hierarchy, which we do not consider in this section; this means that step 3 in Sec. 12.21 will 
not be necessary). For most of this section, we assume real mass matrices, but we will 
demonstrate the dependence on the Dirac-like phases 5 e and 5". In all of the figures shown 
in this section, a "valid" model corresponds to a model generated by our procedure which 
is consistent with current bounds. Therefore, the interpretation of the figures is different 
from so-called "scatter plots" showing the parameter space density for certain assumptions 
on the input variable distributions. Our figures represent discrete predictions for specific 
models, where each point/model can be connected with a specific texture. 

Note that the interpretation of the distributions of observables as predictions has to be 
done with great care, since in nature only one model is actually implemented. In addition, 
one has to use a measure in theory space, which we have done by our discrete choices of 
mixing angles in TJt and U u . However, one can use excluded regions in the parameter space 
(for the given assumptions) for strong conclusions. In addition, it is useful to study the 
model parameter space as a whole in order to obtain hints on how experiments can test this 
parameter space most efficiently (cf. Sec. I4.4p . 

4.1 Distribution of Mixing Angles for Fixed e 

Let us first of all fix e to 0.2, as we have done before. In this case, our best sample contains 
all allowed models for this fixed e, which are 2 468 out of the initial 262 144 models. We show 
the distributions of mixing angles for these models in Fig.[T]in terms of histograms, where the 
gray-shaded regions mark the current 3<r-excluded regions. The vertical lines/arrows mark 
the 3<T exclusion potential of selected future experiments H For sin 2 26*13, the distribution is 
peaking at rather large values, i.e., most models are within the range of Double Chooz or 
the T2K and NOz/A superbeams. T2HK could then exclude most of the rest, and a neutrino 

3 The sensitivities are, for Double Chooz, taken from Ref. [14], for T2HK, taken from Ref. [53], for the 
neutrino factory, taken from Ref. [54] for two baselines, and for SPMIN/SADO (including solar data), taken 
from Refs. [16,55] for a luminosity of 20GWktyr. 
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Figure 1: Distributions of sin 2 26>i3 (left), sin 2 #23 (middle), and sin 2 612 (right), of the models of our best 
sample for e = 0.2 fixed. The bars show the number of selected models per bin, i.e., per specific parameter 
range. The gray-shaded regions mark the current 3cx-excluded regions (c/. Table [1]). 




0.225 0.25 0.275 0.3 0.325 0.35 0.375 
s\n 2 0, 2 



Figure 2: Distributions of sin 2 28 13 (left), sin 2 623 (middle), and sin 2 9 12 (right), of the models of our 
best sample for the optimal e (see main text). 
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Figure 3: Distributions of sin 2 26 13 (left), sin 2 Q23 (middle), and sin 2 12 (right), of the models of our 
best sample averaged over the Dirac phases 5 l and 8 V using uniform distributions (see main text) . 
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factory almost all of the rest. We will discuss the number of allowed models as a function 
of the expected precisions in greater detail later. In fact, the smallest sin 2 2# 13 in all of the 
best sample models is 3.3 • 10~ 5 , i.e., there is not a single model with sin 2 2# 13 = 0. This 
value is not far below the reach of a neutrino factory]^ 

For sin 2 623, the selected models peak around maximal mixing and at sin 2 623 = 0.5 ± e, as 
one may expect by considering e as a perturbation. This means that the value of e may 
determine the deviation from maximal mixing for a large class of models. Interestingly, the 
peaks around 0.3 and 0.7 are already under strong experimental pressure right now, which 
means that an exclusion may come shortly. However, choosing a somewhat smaller e may 
change this argument, as we will discuss in the next subsection. Note that there are gaps 
between the outer peaks and the maximal mixing peak with no models at all, which makes 
these two classes very distinguishable. We have also tested the predictions for 623 > 7r/4 
versus 623 < 7r/4, but we have not found a substantial deviation from a 50:50 distribution. 

As far as sin 2 9u is concerned, the distribution is rather flat in the currently allowed region. 
In fact, the main peak around sin 2 6 l 12 ~ 0.36 is off the current best-fit value (but well 
within the currently allowed range), which means that future precision measurements (such 
as by a large reactor experiment SADO/SPMIN) could exert strong pressure on the models. 
Except from these parameters, we have also tested the predictions for 5. Since we assume 
real matrices in this part, only and ir can be generated (CP conservation). We have not 
found any substantial deviation from an equal distribution. 

It is now interesting to compare our distributions with the literature, where we have chosen 
two specific examples for sin 2 2# 13 predictions. In Ref . [56] , random mixing matrices and 
their predictions for sin 2 # 13 were investigated ("anarchy"). Our distribution for sin 2 2# 13 
follows the general anarchy trend, which is not surprising for a statistical ensemble of models. 
Note, however, that our initial assumptions are qualitatively very different and we have a 
class of discrete models here. For example, we observe a number of discrete possible models 
for very small sin 2 2# 13 . One can also see this excess on a linear scale in sin 2 2# 13 : The more 
or less uniform distribution in sin 2 2# 13 (creating a distribution in log(sin 2 2#i 3 ) peaking at 
large values) has a tendency to small values of sin 2 2# 13 . Another study, which reviews 
existing models in the literature by the model class and their predictions for sin 2 2#i 3 , is 
Ref. [50]. Although our model is motivated by GUTs, it does not at all show the sin 2 2#i 3 
distribution obtained for GUT models in the literature (which peak around sin 2 2# 13 ~ 0.04). 
In particular, we obtain an excess of very large and very small sin 2 2# 13 models. We find this 
result very interesting, because our approach does not imply a particular matter of taste for 
the result, whereas in the literature, many models might have been biased with respect to 
the outcome. 

4.2 Tuning e 

So far, we have assumed e ~ 8c — 0.2 fixed. However, it may well be that a slightly different 
choice for e will shift specific models in our selection range, i.e., that there will be more 
models allowed. From neutrino physics only, one can derive e using specific assumptions. 

4 Extending the allowed values for sfj and Sy by ± arctan 1 j\[2 would produce models with sin 2 29\^ = 0. 
However, this class of models is not included in our initial hypothesis. 
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For example, for a normal hierarchy of neutrino masses using mi : rri2 : = e 2 : e : 1, 
we obtain from the best-fit values and allowed ranges in Tabled] the range 0.15 < e < 0.22 
(3a) with a best-fit value e ~ 0.18. Of course, these values and ranges are based on more 
assumptions that we have used before, and one may find counter-arguments from the quark 
sector (such as using the Cabibbo angle instead). Therefore, in order to test the effect of 
different e's, we vary e in a symmetric range around 0.2 between 0.15 and 0.25. Then we 
choose, for each model, the e which minimizes our selector value. 

First of all, we observe that the number of valid models according to our selector increases 
from 2 468 to 3 316. We show the distribution of these models in Fig. [2J which is similar 
to Fig. [T]and should be compared with that. We find that for sin 2 # 2 3 (middle panel), the 
gaps between sin 2 #23 — 0.5 ± 0.2 and sin 2 623 — 0.5 (maximal mixing) are getting filled 
by some models. This observation supports the hypothesis that the positions of the outer 
peaks are indeed determined by the value of e. For sin 2 612 (right panel), the main peak is 
now at the current best-fit value sin 2 12 — 0.3. This means that the relatively small error 
on sin 2 # 12 exerts pressure on the choice of e. However, this also implies that an adjustment 
of e can, for most models, circumvent the fact that for fixed e most models do not hit the 
currently allowed value of sin 2 #12 for future precisions of sin 2 612 (cf. Fig. [TJ right panel). 
As the last indicator, we observe that the sin 2 29\3 distribution is now even more peaking 
at larger values of sin 2 2# 13 (c/. Fig. [2j left panel), with most models being very close to 
the current bound. However, a small number of models with small sin 2 2# 13 still survives. 
In summary, we find that adjusting e does not change the qualitative conclusions from the 
previous section, but it allows for even more allowed models with a tendency of better 
sin 2 9\2 fits and larger predicted values for sin 2 29\3 . 

4.3 Impact of Non- Vanishing Dirac-Like Phases 

We now discuss the impact of Dirac-like CP phases 5 e and 5 V in Ug and U v , respectively. 
So far, we have assumed that all phases be or tt in Eq. ( I8"a| . and we have obtained the 
CP conserving values and tt for 5 in t/pMNS without significant preference for either one. 
Let us now still fix tpi, (p 2 , <f>i, and 02, in Eq. ([Hat) to their CP conserving values or tt, 
and average over all possible Dirac-like phase values of 8 and 5 U . Note that this implies 
that there will be non-trivial values for both the Dirac and Majorana phases in ?7pmns, 
which can be interpreted for predictions as well. Since these predictions are a bit more 
model-dependent than most of this section, we discuss them in the next section. 

As far as the procedure is concerned, we now vary 5 e and 5 U from to 2tt in 32 equi-distant 
steps (excluding 2tt), i.e., we generate 32 x 32 = 1024 models instead of four before (two 
values and tt times two phases). This implies that we simulate a uniform distribution in 
these phases. Then we choose the allowed models according to our selector, calculate the 
parameter predictions, and, in order to compare with previous results, divide the number of 
valid models in each bin by 1 024/4 = 256 (averaging). Indeed, using this averaging process, 
we find a number of 1 570 valid models (instead of 2 468 for the Dirac-like phases fixed at 
the CP conserving values). Though the order of magnitude is the same, this implies that 
about 64% of all phases are as good as and tt (on average), whereas one may naively 
expect about 50%. 
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We show the effects of the Dirac phase averaging in Fig. [3j which should be compared 
with Fig. [TJ The qualitative result is, again, very similar to Fig. [TJ For sin 2 2# 13 (left 
panel), the distribution becomes smoother because of the averaging process. Similar to 
the e adjustment, larger values of sin 2 2#i 3 become more pronounced. For sin 2 # 23 (middle 
panel), we find some models in the gaps between 0.5 and 0.5 ± 0.2, as one may expect from 
a smearing of a subset of models. In addition, some peaks at around 0.5 ± 0.2 become 
relatively more pronounced, because there is a subset of models which is independent of 5 e 
and 5 V (for 9[ 3 = and/or 9" 3 = 0). For sin 2 6 I 12 , we observe that this smearing fills the gaps 
which are present for fixed e, but the qualitative picture remains unchanged. Therefore, we 
do not find significant deviations from our earlier distributions. However, there is a slight 
excess of models favoring 9 2 3 < it/4 versus 6*23 > tt/4, as it is already obvious from Fig. [3j 
middle panel. As the most interesting part, we now also obtain non-trivial predictions for 
S and the Majorana phases in £/pmns- We discuss these in the next section. 

4.4 Increase of the Experimental Pressure 

Here, we analyze future improving experimental constraints and how they affect the selection 
of models, where we restrict this discussion to e = 0.2 fixed. First of all, we show in 
Fig. H] the distribution of models in the sin 2 2#i 3 -sin 2 #23 (upper left), sin 2 #i 2 -sin 2 #23 (upper 
right), and sin 2 26 l i3-sin 2 9 12 (lower left) planes, where each point corresponds to one or 
more models with a specific parameter combination. The figure shows all the models from 
our best sample based on current experimental bounds, whereas the dark stars mark the 
models remaining after increasing the experimental pressure on a time scale of ten years. 
The corresponding scenario "In ten years" is defined in Table [TJ where we assume that 
#13 will not be discovered. Obviously, the improving sin 2 2#i 3 bound restricts the model 
space very strongly (left panels), where many models close to maximal mixing will survive 
(upper left panel). However, only very few models close to the sin 2 # 12 best-fit value survive 
(lower left panel). In particular, none of the remaining models can reproduce both the 
solar and atmospheric best-fit values closely (upper right panel). Note that the models with 
extremely small sin 2 2#i 3 will be eliminated as well, mainly because of the pressure from 
sin 2 9 12 (cf. lower left panel). Therefore, we expect that the strongest experimental pressure 
on this model space will come from sin 2 2# 13 and sin 2 # 12 , which are two different degrees 
of freedom since they will be driven by different classes of experiments (such as beams 
and short-baseline reactor experiments versus solar and long-baseline reactor experiments). 
Note that this discussion assumes that sin 2 2#i 3 will not be discovered, whereas a sin 2 2#i 3 
discovery would, depending on the new best-fit value, select a different class of models. 

In order to demonstrate the continuous dependence on experimental constraints, we show in 
Fig. [5] the number of valid models as a function of the sin 2 2# 13 bound (left), error on sin 2 # 23 
(middle), and error on sin 2 # 12 (right). The right edges of the plots correspond to the current 
errors, the left edges to exactly known parameters. The light (yellow) curves fix the errors 
for the other two degrees of freedom to their current values, whereas the dark (red) curves 
assume the other errors from the scenario "In ten years" defined in the caption of Fig. [H 
Comparing the light/yellow curves only, the strongest experimental pressure will be exerted 
by sin 2 2# 13 and sin 2 # 12 , which may each independently reduce the number of models to 
one third on a timescale of ten years. This means that sin 2 2# 13 experiments affect the valid 
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Figure 4: Distribution of models for fixed e = 0.2 for the best sample based on current experimental 
bounds (all marks), as well as remaining models for the bounds expected in ten years from now (dark/blue 
stars only); cf. last column of Table[T]for this projection. The distributions are given in the sin 2 26'13-sin 2 6*23 
(upper left), sin 2 6> 12 -sin 2 6*23 (upper right), and sin 2 2# 13 -sin 2 6* 12 (lower left) planes, where the best-fit values 
are marked by lines. Each point corresponds to one or more models predicting these parameter values. 

models as much as solar and potential long-baseline reactor experiments (SPMIN/SADO), 
and already one of these experiment classes already acts as a strong model discriminator. 
The reason is the generic prediction of the value of sin 2 # 12 from the concepts of maximal 
mixing combined with e deviations, whereas maximal mixing is used as an initial hypothesis 
for our models. In order to estimate the combined potential of different experimental degrees 
of freedom, we show as the dark/red curves the number of valid models for smaller errors on 
the parameters not shown (approximately on a time scale of ten years). Obviously, sin 2 2^13 
is a necessary discriminator to exclude all possible models very quickly, whereas no further 
increase in the sin 2 2#i3 bound would make a further model discrimination very hard. 
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Figure 5: Number of valid models as function of the sin 2 2#i3 bound (left), error on sin 2 6*23 (middle), 
and error on sin 2 6*12 (right). The right edges of the plots correspond to the current errors, the left edges to 
exactly known parameters. The light (yellow) curves refer to the current errors for the other two degrees of 
freedom. The dark (red) curves assume the other errors to be as in the scenario "In ten years" (cf. Table [1]). 
The other indicated experiments marked by the vertical lines are referred to in the caption of Fig. [TJ 

5 Distribution of Observables: Dirac and Majorana PMNS Phases 

We now demonstrate that the results in Sec. 14.31 can be used for non-trivial predictions 
of the Dirac and Majorana phases in C/pmns- We have decoupled this discussion from 
the previous sections, because the results in this chapter will be more dependent on the 
assumptions used. However, we believe that already the procedure employed in this section 
warrants some attention, because it opens a new way of model building predictions. The 
actual dependence on the assumptions, which need to be chosen according to the underlying 
theory, deserves further studies. Note that, as in the last section, the "predictions" can be 
used to study the impact of experiments on the parameter space of models, but one has to 
be careful to interpret them as predictions for the actual theory. 

In Sec. I4.3[ we have assumed that (px, tpi, 0i, and 02, in Eq. (IBal) be or 71", but the Dirac-like 
phases 8 l and 5 V be uniformly distributed. We have then generated all possible models for 
different phases and selected our best sample which has been compatible with current data. 
These phases generate non-trivial Dirac (5) and Majorana (0i, $2) phases in the product 
^pmns — Ujlf u , which can be predicted from our best sample. One can now argue if our 
set of assumptions is plausible (such as the set of values used for the mixing angles), if one 
should include variations of the other phases, if the uniform distributions for the Dirac-like 
phases make sense, etc. Therefore, one has to be careful with the interpretation of the results 
in this section. It is, however, the purpose of this section to demonstrate how our procedure 
can provide non-trivial predictions for phases under certain assumptions. 

First of all, we show in Fig. [H] the distribution for 5 (in t/pMNs) f° r our best sample. As 
before, we have generated 32 x 32 models for all possible pairs of Dirac-like phases S e and 6", 
have chosen our best sample, and then have normalized the histograms by 4/1024 = 1/256 
(because we replace four choices of the phases by 1024). In the left plot, we show the 
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Figure 6: Distribution of the Dirac phase 8 for our best sample for all models of that sample (left), only 
models with at least one 6f 3 > or 9" 3 > (middle), and models with both 6*f 3 > and 9" 3 > 0. There are 
no particular assumptions for #13 (in C/pmns)- 

distribution for all models in the best sample. Obviously, there is a strong peak at 5 ~ it, 
which predominantly comes from models with 9\ 3 = 9 V 13 = 0. In this case, Ug and U v are 
real matrices, and the product has to be real as well, which only allows the phases and 
7r for 5. There is no smearing of these models coming from the phase averaging, because 
5 e and 5 U are undefined. This means that the averaging over these phases leaves the peak 
untouched and smears out many of the other models, leading to a relative enhancement 
at the CP conserving choice. The strong preference of tt compared to comes from the 
selector from a non-trivial restriction of the parameter space in all three mixing angles, it 
is not present before the selection process. 

In some sense, the overall preference of CP conservation is connected with our choice of as 
one of the values for the mixing angles 9\ 3 and 9^ 3 , which means that it is model-dependent. 
One can now argue that it does not make sense to consider the cases with 9[ 3 = or 9\ 3 = 
because in this case the Dirac-like phases are not physical and therefore introduce a bias 
(these models have a stronger relative weight because the phases are not varied). Therefore, 
we show in Fig. [6j middle and right, only the models with at least one "physical" (=defined) 
Dirac-like phase (middle) and two physical Dirac-like phases (right) in Ug and U v . Obviously, 
there is a shift in the preference of 5. However, all of the plots in Fig. [6] have in common 
that maximal CP violation is disfavored. The reason for the peaks close to CP conservation 
is the generation of the PMNS Dirac phase in (t/pMNs)e3 as a combination of the Dirac-like 
phases from Ug and U u (times some factors) plus a constant term acting CP conserving. If 
this constant term is larger than the phase- dependent factors, there will be CP conservation 
disturbed by some e-size (or e 2 -size) contribution coming from the mixing angles. If this 
constant term is smaller than the phase-dependent factors, the uniform phase distributions 
will translate into uniform distributions for the PMNS Dirac phase. Note that we have still 
kept all the other phases (pi, <p2, 0i, and 02, fixed at their CP conserving values - which is 
a model dependent assumption (and a constraint from computation power). However, note 
that 0i and 02 only affect the Majorana phases in £/pmns- 
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Figure 7: Distribution of the Majorana phases 0i and <f>2 from our best sample. Each point corresponds 
to one or more models predicting the corresponding phases. 

Our approach does not only predict the PMNS Dirac phase, but also a non-trivial set of 
Majorana phases. We show in Fig. [7] this distribution for our best sample. Obviously, there 
is some clustering close to the CP conserving values, but, in principle, the whole parameter 
space is covered. In addition, note that the parameter space [n, 27r] is redundant for Ov/3/3 
decay, because only the square of the phase enters. 

The predictions for these phases for one specific model become very interesting in combina- 
tion with the mixing angles for that model, because this set of parameters can be used for a 
direct prediction of the mass matrix element for Ou/3/3 decay for Majorana neutrino masses. 
We follow the calculations in Ref. [57] (for a recent review, see also Ref. [58] and references 
therein), where the absolute value of the mass matrix element m ee , which is proportional to 
the rate of the OvfiP decay, is given b}@ 

with m ee = |m« | + \mf} \ e 2l ^~^ + \mfj | e ~ 2i{ ^ +S) (36) 

5 Note the change in notation for the phases: We read off <p\ and <f>2 in a way which is related to a and 
/3 in Ref. [57] by a = 4>% - <f>i, /3 = -<f>i - S. 
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Fi gure 8: Distribution of \m ee \ from our best sample for a normal mass hierarchy and mi as chosen in 
the plot labels. 
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Given the mass hierarchy and using the mass squared differences from Table [TJ one can 
use these equations to compute | function of the lightest neutrino mass m (mi 

for the normal hierarchy, and m 3 for the inverted hierarchy). Each of our models predicts 
the combination (0i, <p 2 , 0i3> ^12> 8) relevant for decay, which results in a particular 

prediction of | function of m. We show this prediction for the normal mass hierarchy 

for two different values of m = mi in Fig.[8j The left plot corresponds to a vanishing mi, the 
right plot to the "chimney", where |m ee | may even vanish. From these distributions one can 
read off what may be obvious: Making the 0^/3/3 decay rate vanish could mean fine-tuning, 
because it requires a specific combination of phases and mixing angles (cf. Fig. 1 in Ref. [57] 
for illustration). Therefore, none of our models predicts a vanishing Ov/3/3 decay rate. 

We show in Fig. [9] the theory allowed regions (left) and our model predictions (right) as 
function of the lightest neutrino mass. The general theory prediction is obtained by varying 
the mixing angles (and phases) within their current 3a allowed ranges. Because of the 
relatively weak sin 2 26 Vi bound the chimney is not explicitly visible in the left plot, and 
a vanishing 0h>(3/3 decay rate is possible for the normal hierarchy. In the right panel of 
Fig. [U we compute for each m the minimum and maximum of |m ee | from all models in the 
best sample. In addition, we show the curves where 99% and 95% of all models are above. 
As the most interesting result, we find that 99% of all normal mass hierarchy models are 
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Figure 9: General theory allowed regions (left) and distribution of our best sample prediction (right) for 
\m ee \ as a function of the lightest neutrino mass m. In the left plot, we show the currently allowed regions for 
|m ee |, where the mixing angles are varied in the current 3cr ranges from Table[T] (computed with the formulas 
given in Ref. [57]). Note that fixing #13 would result in the appearance of the "chimney". In the right plot, 
the bands correspond to all possible models from our best sample. The different intermediate curves limit 
the regions where 99% and 95% of the models can be found above. The vertical line corresponds to the 
choice of Fig. [SI right (for the normal hierarchy). Note that in the right panel, we use both the mixing angle 
and Majorana phase predictions for each model, i.e., we use the combination of all available predictions 
for that model (and not only the Majorana phases). In both panels, the limit from cosmology is obtained 
from a combined SDSS and WMAP analysis [10], and the the limit from Of/3/3-decay is obtained by the 
Heidelberg-Moscow collaboration [59] (with an uncertainty coming from the calculations of the nuclear 
matrix elements). We fix the mass squared differences in both panels to their best-fit values in Table [TJ 



above ~ 0.001 eV independent of ni\. In addition, comparing the two panels of Fig. [9j we 
find that the actually predicted ranges will be much more narrow than the corresponding 
theoretically allowed ranges. Again, this result depends on the assumptions used for the 
phase generation, in particular, the selection of values for the mixing angles, the uniform 
distributions of the Dirac-like phases 8 and 8 U , and the selection of CP conserving values 
for the other phases. However, the obvious observation that one needs to fine-tune the 
parameters in order to make 0z//?/5 decay vanish, seems to be quite general. Therefore, 
we expect that in most concrete realizations, there will be no n- vanishing Oz//5/5 decay for 
Majorana neutrino masses. 
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6 Summary and Conclusions 



As a minimal unified approach to the fermion mass and mixing parameters in the Standard 
Model, it is attractive to assume that all deviations from symmetries or zeros, such as 
deviations from maximal atmospheric neutrino mixing, the small mixing angle 6*13 in {/pmnSj 
the Cabibbo angle 9q in Vckmj an d the mass hierarchies, may be described by powers of 
a single small quantity e ~ 9q. This small quantity can be motivated by Grand Unified 
Theories connecting quarks and leptons. Quark-lepton complementarity can be interpreted 
as a phenomenological implementation of this approach, obtaining the solar mixing angle as 
#12 + #c — 7r /4- We have introduced an extended quark-lepton complementary approach as 
an extension of this relationship. Since £/pmns arises as a product of the charged lepton and 
neutrino mixing matrices, i.e., C/pmns — UjU u , we have postulated that all mixing angles 
in Ui and U v be given by either maximal mixing or powers of e. In this way, the observed 
solar mixing angle can only result from taking the product in £/pmns = U\U V . In addition, 
deviations from maximal atmospheric mixing, # 13 , etc., emerge as predictions in this minimal 
unified approach as being given by powers of e. Note that we have not assumed specific 
forms for Ue and U u , such as bimaximal mixing. However, we have obtained configurations 
involving bimaximal mixing as special cases. 

Using this assumption for the mixing angles in Ue and U„ up to the order e 2 , i.e., the 
mixing angles be 7r/4, e, e 2 , or 0, we have systematically tested all possible CP conserving 
configurations, a total of 262 144 generated models. Naturally, we had to introduce an 
automated selector for the search for models compatible with current experimental data, 
leading to a best sample of 2 468 models. As a first analysis, we have scanned for particularly 
interesting textures, such as with Ue — Vqkm and U v ~ t/bimax, or models being perfectly 
consistent with current data. We have found that these models imply that #13 be large, which 
means that many of these models could be ruled out very soon due to stronger bounds on 
9 13 . However, if 9 13 turns out to be large, we have identified examples with (#12, #13, #23) = 
(33.4°, 7.5°, 43.5°), providing a perfect fit. As the next step, we have particularly focused on 
models with small #13 and increased pressure on the other oscillation parameters in order 
to further reduce our sample. These models are interesting because they could survive 
the next ten years even if 9 13 was not discovered. We have classified the corresponding 
textures systematically by their leading order entries, and we have found 20 distinctive sets 
of textures. These textures can be divided into different classes, such as lopsided models, 
anarchic models, etc. As a very interesting class, we have identified a number of models with 
maximal mixings 9[ 2 , 9 23 , and 9\ 3 , which we have called "diamond models" (because of the 
diamond-like structure of the small entries in M^ aj ). We have shown that the sum rules 
for this and the anarchic-like classes of models are qualitatively different from the ones in 
standard quark-lepton complementarity. 

In an independent approach, we have investigated the mixing angle distributions for our 
best sample compatible with current data. These distributions can be used to study the 
impact of future experiments on the model parameter space, and should be interpreted as 
predictions in that way. As the main result, we have found rather large values of sin 2 26> 13 
preferred, and sin 2 29 13 > 3.3 • 1(T 5 for all of our models. In addition, sin 9 23 peaks around 
0.5 and at 0.5 ± e, whereas sin 2 9 12 peaks above the current best-fit value. Compared to 
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the GUT model literature (surveyed in Ref. [50]), we have not found a characteristic peak 
of models at around sin 2 2# 13 ~ 0.04. Since our matrix generation and selection process 
has been quite un-biased, this could point towards a bias in the building of specific GUT 
models. Compared to anarchic models, we find an excess for small sin 2 26> 13 . In the future, 
especially sin 2 2#i 3 experiments (such as superbeams or short-baseline reactor experiments) 
and potential sin 2 9u experiments (such as a long-baseline reactor experiments) will put 
the model parameter space under pressure. Because we have used maximal mixing as an 
input and have generated deviations from that by powers of e ~ 9 C , sin 2 9 12 is an important 
discriminator for this class of models. However, in order to exclude all models, very strong 
bounds for sin 2 2# i3 are needed, such as those coming from a neutrino factory. 

In a more specific predictability part, we have generated 5 e (in Ug) and 6" (in U v ) with 
uniform distributions in order to obtain predictions for 5 and the Majorana phases in ?7pmns- 
This set of assumptions is more model-dependent than the rest of this work, but it allows 
a very powerful handle on phase predictions. For example, we have found that for 5 (in 
£/pmns), maximal CP violation is significantly disfavored, because it requires large imaginary 
parts that are hard to obtain from the construction £/pmns — U\U V . Furthermore, we 
have combined the non-trivial predictions of the Majorana phases with the mixing angle 
predictions for each model, and we have predicted the Qv/3/3 decay rates. Not a single of 
our models predicts a vanishing Qv(3(3 decay rate, because the necessary phase cancellation 
requires fine-tuning not present in the discussed model parameter space. We have found 
that 99% of all models predict \m ee \ > 0.001 eV independent of the mass hierarchy and 
lightest neutrino mass, i.e., the "chimney" is, in practice, not present. 

In conclusion, we have used a novel approach for studying neutrino mass matrices, which is a 
mixture between basic fundamental assumptions and the systematic machinized parameter 
space scan of a very high dimensional parameter space. This approach has turned out to be 
extremely powerful, because it does not require the diagonalization of matrices. The primary 
objective of this work has been to stay as far away from specific assumptions - which may 
introduce a bias - as currently possible from the computational point of view. With this 
approach, we have not only been able to scan a large parameter space systematically, but also 
been able to make predictions based on a large sample of models compatible with current 
data. Each of these predictions for the mixing angles can be connected to a particular 
texture, which is different from parameter space scans using particular assumptions for the 
input variable distributions and investigating the parameter space density. Only in the last 
part, we have combined these two methods because we have used uniform distributions 
for the Dirac-like phases as an assumption. Naturally, we have not been able to present 
all of our results here, but have focused on the most interesting parts - in our opinion, 
of course. The interested reader can find a tool to view all models from our best sample 
compatible with current data in Ref. [60]. Finally, we believe that the connection between 
quarks and leptons could be the key element in the motivation of future neutrino facilities, 
and we have demonstrated how this element can be implemented phenomenologically in a 
straightforward scheme. 
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